Complexity of the Sherrington-Kirkpatrick Model in the Annealed Approximation 
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A careful critical analysis of the complexity, at the annealed level, of the Sherrington-Kirkpatrick 
model has been performed. The complexity functional is proved to be always invariant under 
the Becchi-Rouet-Stora-Tyutin supersymmetry, disregarding the formulation used to define it. We 
consider two different saddle points of such functional, one satisfying the supersymmetry [A. Cavagna 
et al, J. Phys. A 36 (2003) 1175] and the other one breaking it [A.J. Bray and M.A. Moore, J. Phys. 
C 13 (1980) L469]. We review the previews studies on the subject, linking different perspectives 
and pointing out some inadequacies and even inconsistencies in both solutions. 



INTRODUCTION 



The organization of thermodynamic states in complex systems, fundamental in the understanding of the dynamic 
properties, is a rather difficult task to deal with. The quenched disorder characterizing these systems and the 
consequent frustration are such that a huge number of stable and metastable states arises, even growing exponentially 
with the number of elements composing the system. In order to throw light on the structure of the landscape of the 
thermodynamic potential, a very important theoretical tool is the logarithm of the number of states, either called 
complexity or configurational entropy. 

In structural glasses, seen as disordered amorphous solids and, thus, treated by the techniques of complex systems, 
at the glass transition the entropy drops to a (often much smaller) value going from the liquid to the solid state and 
the states, as opposed, e. g., to the crystal states do not display any specific symmetry. The condensed system has 
lost the ability of visiting different states (at least on the observation time-scale considered) and this codifies into the 
loss of entropy. All the other possible states, not selected at the moment of the transition of the liquid to a glass, are, 
anyway, still there from a statistical point of view and could, in principle, still be reached on much larger time-scales. 
' The complexity counts the many equivalent states that could have been chosen at the moment of the quench. 
. In disordered mean-field models for glasses, e.g. the p-spin interaction spin glass models with stable one-step Replica 
£SJ ' Symmetry Breaking (1RSB) frozen phase [1-3], the presence of many metastable states can be detected looking at the 
relaxation that displays a dynamical transition, with diverging relaxation time scales, at a temperature (the dynamical 
temperature) above the point where a true thermodynamic transition takes place, called static temperature. In the 
range between the static and the dynamic temperature the complexity displayed turns out to be extensive. 

What happens moving from IRSB-stable spin glass models to spin glass mean-field models whose frozen phase is 
| described by means of a Full Replica Symmetry Breaking (FRSB) solution? 

In the last twenty years basically two different proposal have been put forward for the nature and the behavior 
of the complexity of the spin-glass mean-field Sherrington-Kirkpatrick (SK) model [4], the prototype of mean-field 
spin-glass models. The first one was originally presented by Bray and Moore [5]; there the complexity was analyses 
both in the annealed approximation, i.e. as the logarithm of the disordered average of the number of states, and as 
quenched average, i.e. the average of the logarithm, performed making use of the replica trick in the case of replica 
symmetry. The second one was initially proposed in Ref. [6], where Parisi and Potters showed that the complexity 
could be obtained by calculating the partition function of m distinct real replicas of the system [7] and provided the 
connection with the previous formalism by means of the generalization of the two-group Ansatz [8]. In that context 
the annealed solution was called 'unbroken' two-group. 

Over the years it has become more and more evident that an important role in the study of the complexity of 
disordered systems is played by the so-called Becchi-Rouet-Stora-Tyutin (BRST) symmetry. Such a property, first 
discovered in the quantization of gauge theories [9] is a sMpereymmetry, in the sense that it transforms bosons into 
fermions and vice- versa. In the context of stochastic field equations, it can be shown that the integration of the 
generating functional of correlation functions over disorder leads to an action presenting BRST symmetry [10] (for 
the random field Ising model case see also Refs. [11]). The integrated generating functional formally coincides with 
the average over the quenched random couplings of the number of metastable states of mean-field spin-glass models. 
In such a context, the property of BRST-susy invariance has recently been analyses in Refs. [12,13]. In this approach, 
imposing the invariance helps in simplifying the computation and it is equivalent to set the due interdependence 
between the physical objects deriving from the TAP free energy functional and composing the action, namely the 
TAP functional itself, its first derivatives with respect to the average site magnetization (i.e. the TAP equations) 
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and its Hessian. In Ref. [13] a BRST symmetric annealed complexities is built, while the BM complexity functional 
appears to violate such symmetry. 

Generally speaking, the fact that the solutions of a set of equations are BRST-invariant is connected with some 
robustness of the equations under linear perturbation. Let us take into account the equations d mi F = 0, where F is 
the thermodynamic functional and the derivative is taken with respect to the microscopic variable m„ i = 1, . . . , N 
and let us assume that the functional computed over the solutions takes If we perturb the equations by means of small 
external fields, d mi F = — > d mi F = hi, the whole set of solutions could, in principle, drastically change. In general, 
solutions will appear/disappear at different values of F with no given prescription for the relative transformation of 
the complexity function and its dominion. The BRST susy can be recovered assuming that the number of solutions 
at a given value of F does not change In Sec. ID we present an argument the explain such correspondence. 

Both the approach of Ref. [5] and the one of Ref. [13] identify an extensive complexity, computing the number 
of solutions of TAP equations [14] in the annealed approximation and with the further simplification of neglecting 
the modulus of the determinant of the Hessian of the TAP free energy functional. However, the two complexities 
display many differences. They have a different lower band-edge value of the complexity curve vs. free energy and a 
different dominion in free energy, with different thresholds. Their magnitude is different, BM complexity being much 
larger than the one satisfying the BRST-susy. Performing the integration over the whole interval of possible free 
energy values yields a finite result in the BM case, but zero in the BRST-susy one 1 . At the threshold, the behavior 
is once again qualitatively different: the BM complexity goes to zero smoothly, whereas the BRST-susy one drops 
discontinuously to zero. 

What do the two different complexities represent and which one of the two is the "right" one, yielding proper 
information over the organization of the states? Moreover, what do these quantities become in an exact FRSB 
quenched computation? 

We study the complexity of the Sherrington-Kirkpatrick (SK) model, mean-field spin model with quenched disorder, 
critically reviewing the analysis made at the annealed level in the far and near past, linking apparently different 
approaches and discussing the role of BRST supersymmetry in this contest. We will carefully look at the limitations 
of the annealed approaches, in particular from the point of view of physical stability, including the incompleteness of 
the one-step RSB approximation for the SK model. 

Our aim is to present here a comprehensive picture of the annealed level, leaving for elsewhere the study of the 
quenched (FRSB) complexity [15]. 

Before dealing with details and derivations we, now, very briefly anticipate the main results of this paper and we 
outline the scheme of their presentation. 

The BM formalism [5,16] is equivalent to the super symmetric one as presented in Refs. [13,12]. Both at the 
'microscopic' level (site commuting and anti-commuting variables) and at the 'macroscopic' level (order parameters) 
the actions in the two formalisms are related to each other by a simple change of variables. The form of BRST 
transformations and of the Ward identities are also obtained in the BM notation. We, thus, can consider just one 
(supcrsymmetric) functional generating one set of saddle point equations. 

Both the BRST-susy and the BM solution (which we will often refer to as BRST-susy-breaking solution) for the 
annealed complexity are solution of the same set of saddle point equations. 

In Sees. IB, IC we derive such properties, comparing the two approaches in all essential details. 

Afterwards, we give some hints in order to understand the physical meaning of the BRST-susy, showing, in Sec. 
ID, the connection between such property and the property of non- bifurcation of TAP-solutions of linearly perturbed 
TAP-cquations with respect to the unperturbed solutions and the BRST-susy. In Sec. I E the identification between 
the annealed BRST-susy solution and the static 1RSB solution at zero magnetic field, main topic of Ref. [13], is 
rederived in the case of the BM formulation. Moreover, in Sec. IF, we recast everything in the formulation of 
the generalized two-group Ansatz of Ref. [6] and we show that breaking the supersymmetry amounts to consider a 
non-trivial Ansatz in the replica calculation. 

In the second part of the paper, section II, we discuss several issues connected with the problem of selecting the right 
solution between the BM and the BRST supersymmetric one. In particular, we show that the sign of the determinant 
of the Hessian cannot be determined by the saddle point solution at the leading order in N and, applying Plefka's 
criterion to the analysis of the physical meaning of TAP solutions, we explain why the parameter B entering the 
determinant has to be set equal to zero, as in Refs. [5,6,12,13], and why, as opposed to what is stated in Refs. [5,17], 
such a choice does not guarantee the positivity of the determinant of the Hessian of the TAP free energy functional. 

We also recall the results of Kurchan [18] about the spontaneous BRST-susy breaking taking place on the BM 



1 This would be the total number of solutions if the modulus of the determinant were not neglected. 
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saddle point: the BRST-susy is recovered by the analysis of the prefactor of the exponent of NY^bm showing it to be 
zero on the BM saddle point. This result, while ensuring mathematical consistency, poses the problem of which is the 
correct prefactor of the saddle point when the modulus of the determinant is not dropped. 

As a byproduct, dealing with the problem that the true thermodynamic description of the low temperature phase 
of the SK model is yielded by a FRSB solution instead of a 1RSB one, in Sec. II D we show that the complexity 
proposed by Bray, Moore and Young [19], as quenched analogue of the above broadly mentioned annealed complexity 
is computed over a BRST-susy saddle point. 



I. COUNTING THE TAP-STATES: THE COMPLEXITY IN THE ANNEALED APPROXIMATION AND 

THE BRST-SUSY. 



Before discussing the properties of the number of solutions of the mean-field equations for the SK model, we very 
shortly summarize the basic, widely known, features of the Thouless- Anderson-Palmer (TAP) [14] formulation. The 
TAP equations for the local average magnetization nii of the ith spin are 



nrii = tanh (hi — f3{l — q)rrii 



Vi = l,...,JV 



'*3 3 



where q is the self-overlap of TAP configurations: 



III; 



(1) 

(2) 



(3) 



and Jij are distributed as 



The associated free energy functional is 



Fta P ({m}) = £ tap ({m}) - T5 tap ({™}) 

N ^2 



£ta P ({m}) = - Jtjm.mj - — (1 - qf 

{0 1 1 | 1 

lo g2 - — (1 - if - g log ( X ~~ m «) + m ' tanh 1 mi | 



(4) 

(5) 
(6) 

(7) 



where £"tap and Stap are, respectively, the expressions for the internal energy and the entropy of a {mj configuration 
and Eq. (1) is yielded by differentiating Eq. (5) [14]: d mi Ft ap = 0. Furthermore one defines the disorder-independent 
"field" 



5K; q) = i tanh 1 m, + (3(1 - q)m,i 



(8) 



so that Eq. (1) can be reformulated as 



d mi Ft a .p=g{m i ;q)-^2Ji j m j =0 (9) 

3 

Using Eq. (9), the expression for the energy of a TAP solution {m| o1 } can be written without making use of the 

s 

E({m^}) = - -L £ g{mi - q) m, - f (1 - qf = -± £ m, tanh" 1 m, - f (l - q 2 ) (10) 



disorder as 



2/3 
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Combining this last result with the expression (7) for the entropy one is able to formally rewrite the TAP free 
energy as a sum of single site-free energies of TAP solution, as 

f({™ so1 }) = E /K; «) = jf E {- lo s 2 - t (1 ~ q2) + \ [mi tanh_1 mi + log (1 _ m ^ } 

where f(rrii, q) is the site free energy computed on a solution of Eq. (9) 



A. The complexity in the annealed approximation 

What we are interested in, is the number of solutions of the TAP equations at different free energy levels, which 
can provide substantial information in order to understand the free energy landscape and, therefore, the organization 
of the states. If we label each of the N solutions with the index a = 1, . . . ,Af, the number of solutions having free 
energy density equal to / is 

N 

p s (/) = E^[^tap({m«})-iV/] (12) 

a=l 

This can be formally transformed as 

AT i jv 

Ps(f) - E / II dm * 6 K - O 3 [ F ta P ({m}) - N f] (13) 

a=l J" 1 i=l 
,1 N 

= Y[d mi 5(d mi F tap ({m})) \det(d mi d mi F tap ({m})\ 6[F tap ({m}) - Nf] (14) 
J — i a — i 



where the Hessian takes the form 

d mi d mj F ta p({m}) = -Jij 



(31 — m 



hi (y ) d- r )) 



Terms of order 1 /N will be neglected since they are not relevant for the present discussion. 

At this point, one can compute the annealed complexity, i.e. the logarithm of the average over the disorder of the 
density of TAP solutions: 



E = -logp.(/). (16) 

where the overline represents the average over the distribution of the Jjj 's. Details of such computation for the SK 
model can be found both in the original paper of Bray and Moore (Ref. [5]) and in the supersymmetric formulation 
of Ref. [13], as well as in many other works, e.g. Refs. [16,20,21]. Here, we only stress the basic fact that, in 
both procedures, the determinant of the Hessian is taken without absolute value, thus implying that the quantity 
actually computed would coincide with the "true" complexity (even though in the annealed approximation) only if 
the determinant of all the solutions we are counting were positive. Such an assumption is hard to verify in the SK 
model. 

B. BRST-invariance of the density of TAP-solutions. 

The density of TAP solutions can be written making use both of anti-commuting {ip}, {tp} and commuting variables 
as [12,13] 

Ps(f) = J VxVmV^ViP e ^(m,x,^) ^ 
S(m,x,if>,ip) = ^Xtd^Ftapdm}) + E i>i^jd mi d mj F tap ({m}) + u [F tecp ({m}) - Nf] (18) 
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where Va = prefactor x J\ i da,. The transformation under which such an action is invariant is [10,12,13]: 



Srrii = e ipi 



5xi = e u tpi 



Hi 



-e xt 



6^ = 



(19) 



The Ward identities generated by imposing the invariance with respect to the above transformations of observables 
are, e.g., 



(ipitpi) = - (m l x l ) 



(20) 
(21) 



where the average ((. . .)) is computed over the measure e s . 

Bray and Moore [5] used explicitly, in their computation, the TAP equations (9) in order to simplify the procedure. 
In Ref. [13] it was claimed that such a substitution brought to an action no longer BRST invariant. Such an 
inconvenient is, however, only apparent. Indeed, shifting the integration variables {x} the BRST invariant form is 
readily restored. The action yielded by the BM procedure is different from the so called susy one because of the 
insertion of Eq. (9), which means 



S(x, m, ip, ip) = Sbm(x, m, ip, ip) + u 



ljs^2g(mi;q)mi 



2(3 



(22) 



Since the quenched disorder enters the action as {xirrij +ip i tpj +u/2 mirrij), a shift on the integration 
variable 



X% > X^ X{ -\- ~1Tli 



is enough to make both actions coincide: 



S(x, m, ip, ip) = S B m(x + -m, m, ip, ip) 



(23) 



(24) 



In the new set of variables {x',m, ip,ip} the transformation (19) keeping a BSRT-susy function invariant remains 
unchanged. The shift of Xi is performed over an integration variable, thus without affecting the final result, yet the 
transformations given in Eq. (19) and the Ward identities computed with the BM measure in the original set of 
variables are changed. 

In the old set of variables the BRST transformation read (as discussed also in Ref. [22]) 



Sxi = e - tpi 



Srrii = e ipi 

and the averages over the two actions are connected by 

(a{x)) = (a(x- |m) 



Hi = -e Xi 



Sipi = 



BM 



(25) 



(26) 



where the average ((• • -))bm ^ s computed over the measure e 5BM so that Ward identities computed with the BM action 
become 



(^)bm = 

w (^)bm 
Inserting Eq. (27) in Eq. (28) one gets 



{miXi) BM + -q 
( x ?)bm~ u ( to ^*)bm + x 9 



(27) 
(28) 



(^)bm = Tl 



(29) 



Even if in the notation of Ref. [13] it seemed that the action Sbm was not satisfying the BRST relations, this was 
exclusively due to the fact that such relations in the BM notation read differently. 
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Moving to the macroscopic level, where the average number of solution is expressed as a function of the parameters 
q, A, A and B, the two approaches continue to be linked and equivalent to each other. The final expressions for the 
complexity in Refs. [5] (BM) and [13] (CGPM) are: 



BM functional/^ (/; u, B, q, A, A) = ext exp {-A q - u f - (B + A)(l - q) (30) 

{u B q A A} 



(B 2 - A 2 ) 



f 1 dm / 



(Am — tanh 1 m 



2 



IB exp— i — / — + \m 2 + uf(m;q) } (31) 



y f 2^ 2 ~q \l-m 2 ) * { 2p 2 q 

B 2 1 1 

f(m;q) = - log 2 - — (1 - q 2 ) + -mtanh~ 1 m + - log(l - m 2 ) (32) 

CGPM functional^ (/; u, B, q, A, A) = ext exp {-A q - (3 u f - (B + A)(l - q) (33) 

{u B q A A} 



{ B 2_ A 2 } 

2p 2 4 



+ - - —u V - A U q - f3 2 u q{\ - q) + log/ 

. 2 



dm / 1 \ i (A m — tanh 1 m) , , 

' IB exp — 2 >- + \m 2 + u<$> {) {m ]q ) ) (34) 



V / 2^ Vl-™ 2 J 1 \ 2p 2 q 

(3 2 1 
$o(n^; ?) = — log2 — (1 — q) 2 + mtanh -1 m + — log(l — m 2 ) 



(35) 



The first expression is obtained making explicit use of Eq. (10) for the TAP energy, while the second one is 
computed without ever using such a relation. Apparently the two expressions differ. Since, however, both actions 
describe the evolution in the parameter-space of the same system, the two formulations must coincide (exactly as the 
'microscopic' description and the final results do) introducing a proper change of variables. Indeed, to link the two 
approaches the following transformation can be set: 

P 2 

A B m = Acgpm + -j-u q (36) 



2 



Asm = ^cgpm + A C gpm + tt u2 <7 (37) 



Notice that the two functions f(m; q) (in the BM formalism) and &o(m; q) (in the formalism of Ref. [13]) are related 
to each other by 

1 1 P 2 

f(m; q) = $ (m; q) - ^ m tanh -1 m - y <?(1 - q) (38) 

Since the formulations are equivalent, the notation to adopt is no more important. We will use the original notation 
of Ref. [5], not simply because it is the eldest one, but rather to have a more direct comparison with other works, 
e.g. the one of Parisi and Potters [6], that showed that the BM action can be obtained from the Legendre transform 
approach of Ref. [7] making a two-group Ansatz on the matrix Q a b entering the computation of the free energy of the 
coupled replicas (in this way they were able to find new solutions of the BM equations called 'unbroken two-group'), 
or the quenched computation of the complexity of Refs. [23,19]. For simplicity we drop the subscript BM in the 
following. 



C. Saddle point equations 

The variational equations, yielding the saddle point (s) values of the parameters for the annealed complexity take 
the form: 

By 

-g^=0 -/ = (/(m;g)) (39) 

<9£ a n r> a2 ( 1 I l — m 



aB =° - j -T-«- Wb(i-.*) >' (40) 



G 



where the average 



a v i ((Am — tanh 1 m) ) a2 

= * + w? (41) 

gf=0 ^A = -^(l- (7 ) + -(mtanh- 1 m) (42) 

^ = - g = <m 2 ) (43) 



i r 1 

(O(m)) = dm 0(m) e £ (™;«^A,A) (44) 

7 = y dm e £("»;«.9.A,A) ( 45 ) 



is taken over the action 



/ 1 \ 1 , o N (Am - tanh 1 to) 2 „ , % 

£(m; u, <?, A, A) = log f ^— — ^ + Bj - - log (/3 2 2^) - A i- + A m 2 + u /(m; q) (46) 

As noted in Ref. [13] the last term of Eq. (41) was missing in Ref. [5]. The free energy of the TAP solutions f(m;q) 
is the one expressed in Eq. (11) or (32). 

Fixing u, therefore leaving / as a free parameter, the above equations have at least two different solutions 2 . One 
solution satisfies the two relations of BRST supersymmetry as stated in Refs. [13,12], which we rewrite here in the 
present notation as 

2 

B+A=-—uq (47) 

A = (48) 

The saddle point Eq. (40), substituting Eq. (43) in it, admits solution for B > (see Appendix A for details). In 
Sec. II A we recall that a general criterion, formulated by Plefka [31], can be applied as a necessary condition to select 
physically relevant solutions. We anticipate that this criterion requires B = 0, so that Eq. (47) becomes a condition 
for A alone. 

In Figures 1, 2 and 3 we show the behavior, at T = 0.2, of both solutions, both versus / and u. The annealed 
complexity computed over the supersymmetric solution goes to zero smoothly as u — > 0~, it displays a maximum 
at some u max (or / t h, if the behavior vs. free energy is considered) and crosses the w-ax at some uq such that 
uq < "max < 0. Unlike the curve S a (/) is not univocal: it displays a cusp at / t h (see Fig. 1) and then turns 

back. 



2 In practice, when we perform the numerical resolution of Eqs. (41)-(43) we can see that there are two solutions of the saddle 
point equations to which our computer program converge and the selection of one solution instead of the other only depends 
on the initial conditions we give. 
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f f 

FIG. 1. Comparison of complexities yielded respectively by the BRST-susy solution (dashed line) and the 
BRST-susy-breaking solution (continuous line) of the annealed saddle point equations at T — 0.2. The range of free en- 
ergy values of the solution breaking the supersymmetry of the TAP-action (46) is much larger than the /-range of the susy 
one. The maximum complexity for this solution is also much higher than the one for the BRST-invariant solution (the value 
at the cusp of the dashed curve). Such a difference increases by lowering T. 

FIG. 2. Detail of the complexity curves for low free energy values. The asterisks stand for the complexity of the static 
one-step RSB solution at zero magnetic field, coinciding with the BRST-susy saddle point. The values of free energy at which 
the complexity reaches zero (sub-exponential growth of solutions of TAP equations with increasing size) are /bm = —0.7693, 
/susy = —0.7651, both below the true equilibrium value at T = 0.2: / cq = —0.7594 (as computed, e.g., in the FRSB scheme). 
Notice that the branch on the left-hand side of the BRST-susy solution in not shown. 



The other solution of the saddle point equations, which will refer to as BM or BRST-susy-breaking solution, is the 
one introduced in Ref. [5]. Such a solution does not satisfy relations (47-48), thus spontaneously breaking BRST-susy, 
its dominion in u (or /) is broader than the one of the symmetric solution, the maximum value of S a is for u = 
and, at any temperature, it is larger than the correspondent value on the BRST-invariant solution. The behavior of 
the annealed complexity is univocal both in u and in /. 




FIG. 3. Complexity vs. u (conjugated to /) for both saddle points. The BRST-susy complexity displays a maximum below 
zero (corresponding to the cusp in the / dominion (see Fig. 1)). It goes to zero as u — > 0. On the contrary the solution with 
broken susy has a maximum for u — 0. 

FIG. 4. The quantity T 2 - ((1 - m 2 ) 2 ) is plotted vs. u. Following the Plefka's criterion, if such a quantity is positive the 
solution is physically meaningful (see text for details), otherwise it is not. In the case of the BRST solution, it coincides with 
the replicon (stability eigenvalue of the qi fluctuations in the one-step RSB static solution with no external field), which we 
plot with asterisks. From the plot one can see that Plefka's criterion is satisfied only by the BM solution. This means that 
if the annealed approximation in the complexity computation would be a reliable one, and if no further inconsistency would 
arise, the physically relevant saddle point over which it should be computed would be the one breaking the BRST-susy. 
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D. Physical Meaning of BRST-susy in the context of TAP-solutions. 



In this section we will follow a reasonment allowing for an intuitive explanation of the physical meaning of the 
BRST supersymmetry in terms of a particular behavior of the solutions of the TAP equations. 

The objects of interest of the present approach are sums over TAP states, multiplied by some weight function, 
typically of the form e P uF ( m ). We make the hypothesis that the set of solutions over which we perform the sums 
does not qualitatively change adding a small magnetic field to a single spin k. This means that each solution a with 
associated magnetizations {mf } goes continuously to a new solution in presence of the field hk, defining the functions 
{m?(h k )}. 

Given a generic observable depending on the magnetizations, its white average, i.e. the average value that the 
observable g takes over each TAP solution, is given by: 

s a 

where M s is the number of solutions we are summing over. 

We, now, assume that (a) no solution bifurcates and (b) no new solution appears. This can take place either as 
a global property of all TAP states, or as a property of a restricted subset of dominant states, selected according to 
their weights. Under these hypotheses this means that we may write 

[sMfc))Li = ^£ff(<(>0) (50) 

s a 

As a consequence, for any n, the following relation holds 

d n [g(Mh))U _ 1 \^ d n g(mf(h)) 



1 ^ d n g(mf(h)) 



dh n Af s ^ dh n 

a 

We start considering g(m) = TO/ c e /3 " F ( m ), which is the average of the magnetization with the weight e^ uF , used to 
count solutions of a given energy. In general we have: 



dh k 



(xkmi) (52) 



where the average on the R.H.S. is computed with respect to the action S(x, m, tp, tp) [Eq. (18)]. However, under the 
above hypothesis, this must be equal to the average on each solution of the derivative of mfce^ uF ( m ) with respect to 
hk, which is 

We can drop the second term since, by definition, it is zero on the TAP-solutions, and we are left with 

d 2 F V 1 



Thus we have 



(x k mi) = {ipiipk) (55) 



i.e., for k = I, the first Ward identity for BRST-susy systems, Eq. (20). To obtain the second Ward identity, Eq. 
(21), we must apply the same procedure to the second derivative of the quantity g(m) = e ^ uF ( m ), proportional to the 
complexity. From Eq. (18) we have: 



d 2 [ e ^ F M] sol 
dhidhj 



= P(xiXj) (56) 



Under the hypothesis of Eq. (50), this must be equal to the following quantity: 
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1 d2 e /3«F(m a ) ^ 1 / d 2 F(m Q ) + ^2^2 dF dm l dF dm s 

Af s ~f dhidhj J\f s ~ \ dhidhj dmi dhi dm s dhj 

where the second term is zero over any TAP solution and we are left with 

£ ^SP e/3uF( "° = ^> (58) 

a 1 3 

Thus we obtain 

(xiXj) =u(ip i ip j ) (59) 

The previous BRST relations can be obtained under weaker conditions than those considered initially, in particular 
the above derivation still holds if we allow the onset of new solutions, provided that they appear only outside the 
ensemble of dominant solutions, e.g. at the threshold. Indeed, in presence of a weight e P uF ( m ) ) only solutions of 
a certain free energy count, thus supporting our assumption that, in the dominion of interest, the total number of 
states, Af s , does not change. Instead, the condition of no-bifurcation seems unavoidable to obtain the BRST relations 
by means of this argument. 



/3uF(m a ) 



(57) 



E. Identification of BRST solution with static solution. Link with Replica computation. 

In this section we very shortly recall the connection between the static solution at the one-step RSB approxima- 
tion and the BRST-susy solution of S a [13]. We compare results at given values of the breaking parameter m in 
replica formalism and at given values of the variable u, Legendre conjugated of /, in the framework of the annealed 
computation of the complexity. 

The total replica free energy at one-step RSB is: 

f rcp = -^(l-2q 1 )-^[(l-m)q 2 1 +mq 2 Q ] - J- J Vz \og J Vz x Pl (m, z Q , Zl ) (60) 

with 

pi(m, z ,zi) = [2 cosh ((3z 0y /qo + 0ziy/qi - go)] " (61) 
Self consistency equations take the form 

Qo = j Vz Q {t&nh (f3z y/cfo + f3z ly Jqi - q )) 2 (62) 

qi = J T>z (tanh 2 ([3z ^/q^ + fiz\^fq\ - q )) (63) 



4 (li - Qo) + j Vz " lo § J Vzi P!( m ' z °' zi) - (logpi(m, z , zi)) 



(64) 



with 



Leaving m as a free parameter (thus ignoring Eq. (64)), at a given temperature we can define the complexity of 
the system as the Legendre transform of [3m / rep : 

Si(/) - max [(3m f - /3m / rep (m)] (66) 

m 

with conjugated variables (3m and /, 

0m = ^ (68) 
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Eqs. (67) and (68) yield the relation between / and m. Introducing Eq. (67) in the Legendre transformation of 
(3m / rep [Eq. (66)] one can obtain the following relation: 



r) f r °P 
S 1 (/) = /3m 2d/ 



dm 



*(/) 



(69) 

»(/) 



P 2 

—to 



(li-Qo)+ J Vz ° lo S J Pi(™, z o,zi) - (logpi(TO,z ,zi)) 



In zero external magnetic field q n = 0. In this particular case Si represents the entropy of hidden states that we 
would get with the method of the m coupled real replicas [7] and a formal connection can be established between the 
complexity of the BRST-susy saddle point of the annealed computation and the one step RSB replica free energy. 
This was shown in Ref. [13] [Eqs. (33)-(35)] but the formal connection also holds in the original notation of Ref. [5]. 

Indeed, the expression of the annealed complexity obtained simplifying A and A in the logarithm of Eq. (30) by 
means of the BRST relations Eqs. (47), (48) comes out to be 



^a(.f;u,q) = log p 3 (f;u,q, A(u,q), X(u,q)) = (3i 



?«|-| [1 - 2q + (1 + u)q 2 ] + 1- log J Vz[2cosh(3z^q\ "-/} 

(70) 

It can be easily seen that, putting qo = in the replica free energy, Eq. (60), this is bounded to Eq. (70) by 

Z a (f;u,q)=/3u[rP(q,-u)-f] (71) 
provided that one makes the identification 

m = — u q\ = q (72) 

Furthermore, the saddle points of such a function S a coincide with those of / rcp , if the correspondence f Tep (q, —u = 
to*) = / holds at the static value of the breaking parameter to, i.e. to* satisfies Eq. (64). 

The identification in Eq. (72) is evident for the overlap: q\ is the self-overlap of states in the replica 1RSB framework 
and q is the self-overlap of TAP-states, i.e. they are just different representations of the same thermodynamic 
observable. For what concerns the connection between the breaking parameter to and the variable u conjugated with 
the TAP free energy, we can observe that the derivative with respect to / of Eq. (70) gives 

ff = -fi* (73) 

Simply comparing this with Eq. (68) one can identify — (3u as (3m, the slope of the complexity as a function of the 
free energy. As a matter of fact we are just saying that Eq. (70) can be seen as the Legendre transform of (—(3u f re P) 
with conjugated variables — (3u and /. 

Notice that, as a consequence of the link shown in Eqs. (71), (72), the average defined in Eq. (44) coincides with 
the one defined in Eq. (65). Therefore we used the same symbol. 

The one-step RSB static solution is not thermodynamically stable, as it can be shown computing the Hessian of 
the replica thermodynamic potential with respect to the order parameters variations around their self-consistcntly 
derived values, given by Eqs. (62), (63). Indeed, the eigenvalue associated with the overlap-overlap fluctuations, i.e. 
the replicon, is 

A = 1 - /3 2 (cosh- 4 (/3V^i)) (74) 
Using the identity cosh -4 = (1 — tanh 2 ) 2 , the stability condition A > can be written as 

T 2 > l-2 g i + (tanh 4 (/3VgTzi)) (75) 

This is the analogous in the replica formalism of the Plcfka's criterion for the physical relevance of TAP solutions [31] 
(see Sec. II A) 

As we can see from Fig. 4, the replicon of the one-step RSB solution at zero magnetic field, as a function of u — —to, 
is always negative in the dominion where the complexity is positive, even on the 'static' point /o = /(to*). From 
the same figure, though, we observe that the equivalent quantity computed over the BM saddle point of the annealed 
complexity is such that Plcfka's criterion is always satisfied. This is a necessary condition supporting the possibility 
that the BM annealed complexity could represent the number of states of the SK model, yet not a sufficient one, as 
we will see in Sec. II. 
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F. The Two-group Ansatz breaks the BRST-susy 



In Ref. [8] Bray and Moore introduced the two-group Ansatz in order to solve the instability problem of the replica 
symmetric solution of the SK model. In Ref. [7] Monasson showed how the formalism of Lcgendrc transforms can be 
applied to mean-field disordered models through the pinning of real replicas in a configuration space extended to to 
copies of the system. 

Parisi and Potters [6] explained how the BM action can be obtained through the method of Monasson provided 
that the symmetry between real replicas is broken according to a generalized two-group Ansatz. Being n the number 
of replicas introduced to compute the quenched average (following the standard scheme [4,24]) and to the number of 
real copies, they analyzed 



lim - log Z mn = cxt 

n— o n f 



log Ps (f)-mpNf (76) 



where Z mn is the partition function ofnxm copies of the system. In terms of the replicated matrix-parameter Q its 
average is 



Z mn = ext exp { N 
Q 



— ran 
4 



TrQ 2 ) + log E ex P Y E Q^>t 

{er£} V ab,cd 



(77) 



where the indexes a, b = 1, . . . , n while c,d = 1, . . . , to. The four index matrix can be expressed as the composition 
of n 2 sub-matrices Q a b of dimension to x to of the form: 



Qab = ( 

The matrices Q^ b are further parameterized as 



■ib 



(78) 



Q± b = Q ab ±^ + ^ (79, 



Furthermore is Q a c a = 0. In Ref. [6] the last term was B a b/y 2 . We write C a b — 2-B a h both in order to avoid confusion 
with the parameter B in the S a expression of Sec. I and to obtain more symmetric expressions in following works. 

Eq. (76) is, then, computed making use of the above Ansatz, getting a complexity that can be formally connected 
to the one of Bray, Moore and Young, Ref. [19] through a given change of variables. We will analyze it more explicitly 
elsewhere [15], where we will use such transformation. 

For the time being we are mostly interested in the annealed case, that can be obtained putting in matrix (78) all 
Qab with a b equal to zero. The Ansatz reduces, then, to consider n diagonal blocks of to x to matrices (78), for 
b = a, where the elements are built with: Q aa — Q, A aa = A, C aa = C. In this case, the change of variables we need 
to connect the BM formulation of the annealed complexity to the Parisi-Potters (PP) one, are 3 : 

BM PP 

q = Q (80) 
A=[A+^Q) (81) 

A = ^ (c + mA+^o) (82) 



2 V 4 

U = —TO 

Writing the equations with the substitutions (80)-(82), allows for an immediate connection between breaking the 
matrices structure into two groups and breaking the BRST-susy. Indeed, Eqs. (81) and (82) transform into the BRST 
relations (47) and (48) if we set A = C = 0, i.e. if we do not break the matrix structure at all (unbroken two-group). 
On the contrary, setting values of A and C different from zero amounts to break the BRST symmetry and leads to 
values of q, A and A independent among each other, i.e. to the BM solution. 



3 Notice that also in Ref. [23,19] there is a Q. In Ref. [19] this is related to q (TAP self-overlap) as Q — [3 2 q. The parameters 
Q in [23,19] (Qbm) and in [6] (Qpp) are proportional to each other but different: Qbm = f3 2 Qpp. 
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II. PROPERTIES OF THE ANNEALED COMPLEXITIES 



In the previous part of the paper, we have studied the so called BM complexity [5] (the one breaking the super- 
symmetry) and the BRST-symmetric complexity [13] as two different solutions of the same saddle point equations, 
derived from the same BRST-susy functional (see Sec. IB). 

In this section we will discuss several issues connected to the problem of selecting a proper, physically meaningful, 
complexity and to see whether any of the two candidates fulfill the minimal requirements. 

We start noticing that not all TAP solutions can be associated to stable thermodynamic states. They have first 
to satisfy the Plefka's criterion (see below), which guarantees the right expression for the linear susceptibility of the 
states under consideration, and, then, to be minima of the TAP free energy functional. 

Operationally speaking, what happens is that, even if on physical grounds we are willing to count only the minima 
of the TAP free energy, in practice, we also count other kind of solutions of the TAP equations. As a consequence, 
it may happen that we have a mathematically correct computation which, however, has no physical meaning at all 
or a rather obscure one. In subsection II A we will see that, at the present stage, we have no way to state that any 
of the two solutions of the complexity saddle point equations correspond to counting only minima of the TAP free 
energy. On the other hand, though it is just a necessary condition, a partial discrimination can be done on the basis 
of Plefka's criterion, which is satisfied by the BRST-susy-breaking solution but not by the BRST-susy solution. 

Since the BRST solution can be linked to the saddle point solution obtained in a 1RSB replica computation in 
zero magnetic field, it is natural to check its stability according to the usual criteria in that framework, e.g. the 
positivity of the replicon eigenvalue; however, as we anticipated in Sec. I E, this condition is formally identical to 
the Plefka criterion, which is, therefore violated by the BRST-susy solution. This violation leads to a mathematical 
inconsistency, as we show in Sec. II A. 

The BRST supersymmctry is a property of the action and, therefore, its violation poses a consistency problem to 
the BM saddle point, too, due to the lack of control on what we are counting. The only possibility for the BRST- 
susy-breaking saddle point to be valid is to guess that it counts only the minima, as it is claimed in Ref. [17]. It 
would describe only a physically stable portion of the set of TAP solution and we would not expect that it satisfies 
global properties valid on the whole set, such as the Ward identities derived from the BRST susy [see Eqs. (20)- (21)]. 
However, this assumption is hard to justify, as we will discuss in the subsections IIB-IIC, recalling the properties of 
the Hessian of the TAP free energy as derived in Ref. [25]. 

The effects of neglecting the sign of the determinant of the TAP solutions, which is at the origin of the BRST 
susy, are particularly dramatic when we set u = 0. In this case we are summing over each extremum of the TAP 
functional with a weight given by the sign of its determinant. The Morse theorem tells us that this quantity is a 
topological invariant and it is equal to 1 in this specific case. While the BRST-susy solution does satisfy the theorem, 
the BRST-susy-breaking saddle point yields a quantity of order e' VEBM 3> 1, if Sbm > 0. 

A first way to recover mathematical consistency in the u = case is to guess that BRST susy is restored, e.g. 
because the prefactor of the exponential is e~ Nj:BM . In this case the expansion of the prefactor in powers of 1/N 
will be zero at all orders (such as the expansion of e -1 / x in powers of x). This series was considered by Kurchan in 
Ref. [18] were it is shown that its coefficients are indeed all null. This proves that the BM saddle point is possibly 
mathematically consistent but rises another difficulty. We already stressed that at the order N it is not possible to 
prove whether the BM complexity is counting only minima or not. When computing the corrections without modulus 
of the determinant of the Hessian one seems to have a vanishing prefactor, though. Indeed, at this stage, in order 
to accept the BM prediction, we must, then, assume that the effect of taking the modulus of the determinant into 
account is to change the value of the prefactor from exponentially small to finite, a completely unjustified assumption 
at the present state. 



A. Plefka's criterion 



In this subsection we recall the results of Ref. [25] on the properties of the Hessian of the TAP free energy on a 
generic configuration {m t }, which is also the inverse of the susceptibility matrix: 

The problem is studied considering the resolvent of the Hessian. As a byproduct of the computation one obtains the 
magnetic susceptibility of a solution. On physical grounds this must be equal to (3(1 — q), however it turns out that 
not all TAP configurations satisfy this relation. Instead the condition for a TAP solution {rm} to yield the right 
physical susceptibility is 
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i 

The TAP solutions which do not satisfy the previous relation turn out to have the following non-physical susceptibility: 

Xi = /?(l-?) + y (85) 

where p is defined as 

P = /3 3 ^£(l-m 2 ) 3 (86) 

i 

Therefore, after having computed the complexity, we must check that the Plefka relation is verified to be sure that 
we are counting physical solutions. 

For N — ► oo, the site-average becomes, in the notation of the annealed computation of Sec. I C, the average of Eq. 
(84) over the action C [Eq. (46)] and the above inequality reads: 

1 - /3 2 (l - 2g + (to 4 )) >0 (87) 

As we can see from Fig. 4 such a condition is satisfied by the saddle point breaking the BRST-susy, but is violated 
by the BRST-invariant one. In the replica language Eq. (87) is the replicon. The identification of the measures over 
which the average is performed is given by Eqs. (70)-(72). 

As we show in Appendix B, the computation of the determinant of the Hessian, which is a crucial step in the 
derivation of the action, is deeply connected to the Plefka computation of the resolvent [25]. As a consequence, one 
sees that, in order to set the parameter B of the determinant equal to zero, as it was chosen in both solutions, one 
must check that Plefka criterion is satisfied. Therefore, the fact that the BRST-susy solution violates it, turns out 
to be not only a problem of physical meaning, but also a problem both for the replica computation and for the 
mathematical consistency of the solution. 

We conclude stressing that although the annealed computation is well defined on mathematical grounds, what really 
matters for physics is the quenched computation where this problem can possibly be cured considering a full-RSB 
BRST-susy solution. As we will see, this certainly happens at the lower band-edge of the quenched complexity (the 
free energy value at which the complexity vanishes) which, as expected, coincides with the equilibrium free energy 
given by the Parisi solution and it is BRST symmetric. 

B. The spectrum of the Hessian of F tap ({m}) 

The spectrum of the eigenvalues of the inverse susceptibility matrix X ( m ) on a generic configuration {m^} is 
determined in Ref. [25] for small eigenvalues and in the region of {to} configurations such that xp ~ 0. The eigenvalues 
distribution is written, at the leading order in N, as 4 : 



p(A)= ^v A ~i' M<<1, A<<L (88) 

The minimum of the spectrum is therefore positive, irrespective of the value of xp, implying that the TAP free energy 
is semi-convex everywhere at the extensive level. In the same paper this result was proven rigorously also for generic 
configurations with a finite xp, either positive or negative. We remark that this result, however, does not exclude the 
existence of solutions with a sub-extensive number of negative eigenvalues. Actually, Morse theorem implies that these 
solutions do exist and, indeed, there must be an equal number of solutions with positive and negative determinant. 
Plcfka's result just states that we cannot distinguish them at an extensive level. 



'in Ref. [25] the factor 1/4 in front of the second term in the square root was missing. 
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C. The determinant of the Hessian 



As we have previously discussed, the fact that no negative eigenvalue exists in an extensive quantity does not mean 
that saddles and maxima do not exist, which is indeed impossible on topological grounds. Therefore, the condition 
xp > does not bypass the problem of the sign of the determinant of x 1 - 

In appendix A the determinant is computed making use of Grassmann variables and using the saddle point method. 
If one evaluates the saddle point over B before evaluating the integrals over the nii, it turns out that there are two 
solutions B = and B > (B is the parameter originally used in Ref. [5] 5 ). If xp > the solution B = satisfies 
the condition for the saddle point of being a maximum of the integrand on the integration path and is the correct 
one, while if xp < the correct solution is the one with B > 0. 

Both the BRST-symmetric and the non-BRST-symmetric solutions assume that B = 0, consequently the expression 
of the determinant of the Hessian of the TAP free energy becomes 

- n {jih?) exp - qf ) m 

Since the prefactor is positive for any {rrii} configuration, this choice of B would provide an a posteriori justification 
for neglecting the absolute value of the determinant if the calculation were exact at all orders in N. Anyway, we 
cannot neglect the fact that all computations are performed at the leading term in N , as N — > oo, thus ignoring all 
sub-extensive contributions. 

As a matter of fact, the sign also depends on neglected prefactors: 



det x _1 = a(N) det x" 1 (91) 

sp 

where a(N) cannot be determined at the leading order in TV, but needs to be obtained from corrections of 0(1): 



1 1 

— log det x- 1 = — log det x _1 



+ lloga(A0 (92) 

sp 1 V 



While the magnitude of the determinant is not qualitatively changed by 0(1) corrections (provided they arc not 
zero, of course), its sign can well be determined by eigenvalues that are present in a sub-extensive quantity. 



D. The Bray-Moore- Young FRSB quenched complexity is BRST-invariant. 

Both the BRST-susy and the BRST-susy-breaking solution give a lower band-edge (the free energy value at which 
the complexity vanishes) different from the equilibrium free energy of the SK model. This is not surprising since we 
are performing annealed averages, while we expect that the physically correct computation is quenched. 

In order to cure this deficiency it is rather obvious that one has to go on computing the quenched complexity in the 
Full RSB scheme, for which the SK model is known to be thermodynamically stable. The quenched case is formally 
discussed in Ref. [22], and one instance of FRSB quenched complexity is presented in Ref. [15]. The first step in this 
direction, however, was performed by Bray, Moore and Young in Ref. [19] where the solution there analyzed had a 
lower band-edge coinciding with the FRSB static one. Now, we notice that the assumptions under which they look 
for such a solution, are exactly the BRST relations (47), (48) generalized to the function order parameters one has to 
introduce in the quenched FRSB scheme of computation. Since we are using their very same notation we can directly 
rewrite here Eq. (19) of Ref. [19], which we break into two lines: 



-P(l) 



A 



1 2 

-rru q EA 



(93) 



5 There, the quenched average determinant of the Hessian was computed as 



det y _1 = ext exp < N 
B I 



B 2 



2 



B(i- q ) + ^(i- q y + -> ;io g ( — n -, + b 



2j3 2 i/ ■ 2 i, ■ N 



^ log (ri 



(89) 



where the initially integer positive number of replicas was sent to —2 at the end of the computation. 
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for the diagonal part (the one surviving in the annealed case), and 




u q{x) 



P 2 2 



u 2 q{x) 



Vx € [0, 1] 



(94) 



for the off-diagonal elements. The further assumption was, then, made 



<?(!) = QEA 



(95) 



The relations expressed in Eq. (93) [inserting assumption (95)] are exactly our Eqs. (47), (48). The same holds for 
the off-diagonal terms if we recognize that, in Eq. (94), the off-equilibrium analogue of A is — p(x) and the analogue 
of A is p*(x)/2. What was found is, then, the quenched improvement of the annealed BRST-susy solution and not 
the quenched analogue of the BM annealed solution of Rcf. [5], for which the above mentioned relations do not hold. 
Moreover, the choice of identifying the elements q(l) of the diagonal block of the FRSB matrix with the elements 
on the diagonal, qEA, leads precisely to the self-consistent, stable, Parisi solution of the SK model [26]. In such a 
case, though, no parameter is left free to vary and, therefore, no analysis over the number of states at given / can be 
performed, i.e. no complexity can be built. 

It is the subject of another paper [15] to go beyond this point and look for a generalized solution that allows a 
'quenched probe' in a free parameter (the state free energy / or the generalization of the m parameter in the Legendrc 
transform approach), yet recovering the right equilibrium value as lower band-edge of the complexity. 



In this paper we have shown that the BM action [5] and the action considered by Cavagna et al. [13] are completely 
equivalent and one can move from one to the other through a simple change of variables. This is also true at the 
microscopic level thus implying that the BM action too is BRST-susy invariant. In particular, this equivalence implies 
that each solution of the BM saddle point equations is a solution of the equations of Ref. [13] as well and vice-versa. 
As a consequence we are left with the problem of selecting the "right" solution. 

In order to determine which one, if any, of the two proposals was the one actually representing the multiplicity 
of stable and metastable states of the Sherrington-Kirkpatrick mean-field spin-glass model we have been critically 
reviewing the properties of both, sometimes solving apparent incongruencies, other times pointing out substantial 
inadequacies. 

We first summarize the case of the BM BRST-susy breaking annealed complexity. The BM saddle point is not 
BRST-susy at any value of u, while we have shown in Sec. I that the BM action satisfies such a supersymmetry, 
although this is differently expressed in the BM notation with respect to the one of Rcf. [13]. As already noticed by 
Kurchan in Ref. [18], we have to be careful and we need some justification before adopting it, if we want to preserve 
mathematical consistency. In particular, one may show that the BRST-susy is restored considering subextensive 
corrections to the saddle point and notice that they amount to a prefactor of exp NT,bm- In the case u — 0, it has 
been shown by Kurchan by means of a series expansion in power of 1 /N that the prefactor of the non-BRST saddle 
point is zero at all orders of the expansion. This could imply a zero prefactor or it could leave the way open to a non 
zero but exponentially small prefactor. Both outcomes, however, strongly change the BM prediction. Therefore, to 
save this prediction, one should show that the effect of keeping the modulus of the determinant amounts to change the 
prefactor from an exponentially small value to a finite one. A step in this direction could be possibly done generalizing 
the technique of Ref. [27] for one dimensional random systems. 

Very recently, in Ref. [17], it has been claimed that the BM saddle point counts only minima. This would imply that 
the complexity does not involve a sum over all solutions and, therefore, is not constrained to satisfy global relations, 
like those imposed by BRST-susy or by the Morse theorem. However, there is no prove at all that the BM saddle point 
counts only minima, since Plefka's analysis of the Hessian shows that all solutions of the TAP equations have strictly 
positive eigenvalues only at the leading order in N. In order to get information about the sign of the determinant one 
should be able to select the minima among all states satisfying Plefka's criterion. This can be explicitly seen in the 
spherical p-spin models where these solutions can be classified as minima and saddles, differing only for one negative 
eigenvalue. The Plefka criterion, xp = 1 — /3 2 (1 — 2q + l/NJ2i m T) > is, thus, not related to the fact that a given 
solution is a minimum or a maximum or a saddle but rather it guarantees that it yields the correct susceptibility. 

The lower band-edge of the quenched complexity, computed with the Parisi Ansatz by Bray, Moore and Young [19], 
gives the correct equilibrium free energy but such solution is not a modification of the BRST-susy breaking annealed 
solution of Bray and Moore [5]. On the contrary, it turns out to be BRST symmetric, as we have shown in Sec. II D. 
This also means that, up to now, no quenched extension of the BM annealed saddle point has been taken into account. 
We will show elsewhere what are the minimal assumptions to break the BRST-susy in the quenched case [15]. 



III. CONCLUSIONS 
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Furthermore, the non-BRST solution is not in agreement with the numerical results recently obtained by Plcfka 
through some modified TAP equations [28] . By means of its method he obtains all the minima of the standard TAP 
free energy satisfying xp > 0. The special set of solutions he collects, i.e. the minima with xp > 0, is precisely the 
set that is supposed to be counted by the BM complexity. Now, the BM solution predicts that the great majority 
of solutions has a value of x p , e.g. at temperature T = .5 given by x p — 0.132408 > (Sbm = 0.002775) and at 
T = .2, by x p — 0.5952975 > (Sbm = 0.05219). Plefka, however, on the basis of its numerical data, hints that all 
the minima of the TAP functional in the thermodynamic limit have a zero x p as N — > oo. Moreover, the number of 
solutions found by Plefka varying N are, e.g. at T = 0.2, less then expTVEeM an d the discrepancy increases with 
increasing N. This is also true for what concerns the range of free energy values over which the complexity is non zero. 
This finding is also in agreement with what has been found at the FRSB level of computation, as will be discussed in 



Notice also that the BM result for xp cannot be changed considering the quenched average instead of the annealed 
one. In Ref. [5] also the replica symmetric quenched complexity was considered and there the authors showed that 
the annealed and the quenched non-BRST saddle points coincide at u > u c (u c < 0) and, in particular, at u = to 
which the highest number of solutions would correspond if the modulus of the determinant of the Hessian was taken 
into account. 

It is worth mentioning, anyway, that the zero temperature limit of the BM total complexity coincides with the 
computation of the number of solutions of the zero temperature limit of the TAP equations, m; = sign(^V Jijirij), 
where no reaction term is present [5,20,21]. If the coincidence of the exactly zero temperature behavior with the T — > 
behavior would be a necessary condition, this would be a strong hint that the BM saddle point provides, indeed, the 
right complexity. However, we notice that this does not exclude the existence of other solutions displaying aT^O 
limit of the complexity different from the value directly computed at T = 0. For instance, in the p-spin spherical 
model, a whole branch of TAP solutions existing at zero temperature disappears as soon as we infinitesimally heat 
the system [29] (see also appendix of [30]). 

Looking at the other solution, we observed that the BRST-susy saddle point does not yield a proper result either. 
Indeed, it counts TAP solutions that do not satisfy the Plefka criterion, i.e. solutions not corresponding to physical 
states. Furthermore, even its mathematical consistency is doubtful, since it is obtained setting B = (see appendices 
A and B), and, according to Plefka's analysis of the resolvent [31,25], this assumption is only justified if xp > 0, 
while for xp < a solution with B ^ must be considered. It may happen that, upon passing to the quenched 
computation, the violation of the Plefka criterion of the BRST solution may be cured. Certainly this happens at the 
Full-RSB lower band-edge, which, as we said above, is BRST supersymmetric [19]. 

In order to solve the problem of selecting a meaningful complexity one can still try to generalize the proof of Kurchan 
[18] to arbitrary u, possibly using a procedure not involving a series expansion in 1/N and/or look for a quenched 
FRSB solution breaking the BRST-susy, thus solving a Parisi equation like the one of coupled replicas investigated in 
Ref. [32] but with a different boundary condition. Both approaches are currently under investigation, together with 
the study of the quenched BRST solution [15]. 

Acknowledgments: We thank A. Cavagna, I. Giardina and A. Annibale for many useful discussions. 

APPENDIX A: THE DETERMINANT OF THE HESSIAN OF THE TAP FREE ENERGY. 

The inverse susceptibility matrix is 



[15]. 



Xij — Jij + O-i^ij 



(96) 




(97) 



where the Jjj are distributed according to 




(98) 



As it is usually done, we will not consider terms of order 1/N in the following computation. 
The determinant can be written with the help of Grassmann variables, (77, 77), as 
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det x 1 = / II dr,i Ar >i exp \ E ^(X-ij^i \ = / II dr,t d?? * exp \ ~ E J v (ViVj + VjVi) + E "^i a ^ 

i—1 I ij J i— 1 | i<j i 



(99) 



Its average over the disordered interaction is 



detx _1 = J Yldrji d^exp j-^ (^ViVi^j +Yl^ a ^ 1 

'J f[ d Vi dfji exp | Vi ( iw + a i) 



(100) 



-w 2 N/2 



-f 

J — c 



du> 
dw 



exp 



exp < ^ log (iw + ai) 



exp (NH{w)) 



H(w)i 



w A 1 



log(iw + a^) 



(101) 



2 AT 

To compute the integral for large N we make use of the saddle point approximation, thus evaluating the solution of 

i 



dH _ J_ V- 

c)w N ^— ' it 



dw " ' N ^— ' + 

d 2 H 
dw 2 



= 



= -i + 4E — 



(102) 
(103) 



The second inequality is a condition that in most cases ensures that the integration path can be modified in order 
to cross the saddle point in the proper way, and is analogous to the maximum condition of the Laplace method. 
Changing iw in 



iv = iw + [3(1 — q) = iw + ^ — 



1 1 



(31 -w 2 



the stationarity condition for the saddle point reads 







N 2-? 1 



(l-m?) J 



and the condition for H becomes: 



+ iv/3 (1 — mj) 



(1-m?) 5 



= 



N^ [1 + ivl3 (i-m 2 )} 

The saddle point equation (105) has two solutions: v = 0, v = v* ^ 

• Solution v = 0. 
In this case Eq. (106) simplifies in 



> 



(104) 



(105) 



(106) 



(107) 



This is exactly Plefka's criterion characterizing a physically relevant solution. The stationary value v — is a 
maximum of the exponent H and corresponds to TAP solutions yielding the physical expression of the linear 
susceptibility. 
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• Solution v = v* . 
For such a saddle point Eq. (106) can be written as 

v2 



^ (i-m?r f3 ^ i-m?r , x 

V . * y „. V i ^ 7 >0 108 

of which the first two terms cancel each other. Since m? < 1 always, this implies that in order to have 



i 3 



*4E (l " mfr 2 >0 d09) 
nZ ~ j [1 + iv/3 (1 - m 2 )f 



must be iw real and positive. In order to go back to BM notation we define the real variable B = iv(3. If B > 
the stationary point B = iv*(3 is a maximum of H. Inserting such a positive value into the saddle point Eq. 
(105) one gets the inequality 

/3 2 (l-m|f 2 ^2 nin . 

z v 1 ' i 

Thus violating Plcfka's criterion. 

Summarizing, if the Plefka criterion is satisfied the correct solution is B = while if it is not satisfied one must 
choose the solution with B > 0. Therefore if one sets a priori B = 0, then he must self-consistently check that the 
Plefka criterion is verified. Therefore, from this point of view, the BRST solution is mathematically inconsistent. The 
Plefka criterion arises as a condition to determine which is the correct solution of the resolvent equation; the fact that 
it is also the condition for determine the correct saddle point for B, as we derived above, is not surprising, indeed the 
two computation are intimately related as we shown in the following appendix. 



APPENDIX B: IDENTIFICATION OF detx 1 SADDLE POINT EQ. AND RESOLVENT EQ. FOR x 1 

In this appendix we would like to stress the analogy between the saddle point equation for w = —i[B/f3 + f3(l — q)] 
(see App. A) and the equation for the resolvent of the inverse susceptibility. Using the notation of Ref. [25], the 
resolvent is 

1111 
R(z) = _Tr r = -Tr (111) 



where the resolvent equation is 



R(z) = ±-Y -!- (112) 

w N ^ z-R(z)-a,i v ; 



and the condition 3(i?(z)) > must hold for ^s(z) < 0. 

We notice that equation (112) evaluated in z = is identical to Eq. (102) in App. A provided the transformation 
5i(-R(0)) = iw is performed. Conversely the condition SJ(i?(z)) > for 3(z) < is equivalent to (103). Indeed the 
derivative of the resolvent comes out to be 

dR - n*) (113) 

^>4E f ~^7T~ ,2 (H4) 
iV (-z + R{z) + a l ) 

The function Y(z) is always positive, for any real z. 

If we expand R(z) around a given real value of z — zr for a small negative imaginary part — it we get 



R(z R - it) = R(z R ) - te Y(z R ) 1 - — 



(115) 
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Thus, the condition on the imaginary part of R(z) for negative 9z (Pastur theorem [25,33]) reads 
lim SSRfz) = e lim Y(z) ( 1 - ^ | > 

Q-z^O- Qz-+0- \ dz J 

Eq. (113) evaluated at z = satisfies condition (116) if Y(0) < 1. 

If we set R(0) — iw the equations (112) and (116) evaluated at z = are equivalent respectively to Eqs. (102) and 
(103), thus legitimating this last equation as a validity condition for the saddle point of H(w). 

The resolvent equation has two roots, which, for small z and in the region where f3 2 (l-2q+ (to 4 )) ~ 1, were; 
evaluated e.g. in Ref. [25]. They correspond to the B — and B > solutions of the previous section. The condition 
Q(R(z)) > for 3(z) < selects one or the other solution depending on the value of xp inasmuch as the condition 
(103) selects the correct solution B = or B > depending on the value of xp 



< 1 



3z=0 



(116) 
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